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Abstract 

The aim of this paper is to study the elastic stress and strain fields of dislocations and 
disclinations in the framework of Mindlin's gradient elasticity. We consider simple but 
rigorous versions of Mindlin's first gradient elasticity with one material length (gradient 
coefficient). Using the stress function method, we find modified stress functions for all 
six types of Volterra defects (dislocations and disclinations) situated in an isotropic and 
infinitely extended medium. By means of these stress functions, we obtain exact analytical 
solutions for the stress and strain fields of dislocations and disclinations. An advantage 
of these solutions for the elastic strain and stress is that they have no singularities at the 
defect line. They are finite and have maxima or minima in the defect core region. The 
stresses and strains are either zero or have a finite maximum value at the defect line. The 
maximum value of stresses may serve as a measure of the critical stress level when fracture 
and failure may occur. Thus, both the stress and elastic strain singularities are removed 
in such a simple gradient theory. In addition, we give the relation to the nonlocal stresses 
in Eringen's nonlocal elasticity for the nonsingular stresses. 

Keywords: gradient theory, dislocations, disclinations, nonlocal elasticity, hyperstress 



1 Introduction 

The traditional methods of classical elasticity break down at small distances from crystal defects 
and lead to singularities. This is unfortunate since the defect core is an important region in 
the theory of defects. Moreover, such singularities are unphysical and an improved model of 
defects should eliminate them. In addition, classical elasticity is a scale-free continuum theory 
in which no characteristic length appears. Therefore, the classical elasticity cannot explain the 
phenomena near defects and at atomic scale. 

An extension of the classical elasticity is the so-called strain gradient elasticity. The physical 
motivation to introduce gradient theories was originally given by Kroner [1, 2] in the early sixties. 
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The strain gradient theories extend the classical elasticity with additional strain gradient terms. 
Due to the gradients, they must contain additional material constants with the dimension of a 
length, and hyperstresses appear. The hyperstress tensor is a higher order stress tensor given in 
terms of strain gradients. In particular, the isotropic, higher-order gradient, linear elasticity was 
developed essentially by Mindlin [3-5], Green and Rivlin [6,7] and Toupin [8] (see also [9-11]). 
Strain gradient theories contain strain gradient terms and no rotation vector and no proper 
couple-stresses appear. In this way, they are different from theories with couple-stresses and 
Cosserat theory (micropolar elasticity). Only hyperstresses such as double or triple stresses 
appear in strain gradient theories. Double stresses correspond to a force dipole and triple 
stresses belong to a force quadrupolc. The next order would corresponds to a force octupole. 

Stress and Hyperstress are physical quantities in a 3-dimensional continuum mechanics. 
Within a framework of the 4-dimcnsional spacetime continuum, stress and hyperstress translate 
into momentum and hypermomentum [12, 13]. 

In the present work, we consider two simple but straightforward versions of a first strain 
gradient theory. We investigate screw and edge dislocations and wedge and twist disclinations 
in the framework of incompatible strain gradient elasticity. We apply the "modified" stress 
function method to these types of straight dislocations and disclinations. Using this method, we 
derive exact analytical solutions for the stress and strain fields demonstrating the elimination of 
"classical" singularities from the elastic field at the dislocation and disclination line. Therefore, 
stresses and strains are finite within this gradient theory. We obtain "modified" stress functions 
for all types of straight dislocations and disclinations. In addition, we justify that these solutions 
are solutions in this special version of Mindlin's first gradient elasticity. We also give the relation 
of nonsingular stresses of dislocations and disclinations to Eringen's nonlocal elasticity theory. 
We show that these stresses correspond to the "nonlocal" stresses. 



2 Governing equations 

Following Mindlin [3-5] (see also [10]), we start with the strain energy in gradient elasticity of 
an isotropic material. We only consider first gradients of the elastic strain in this paper. In the 
small strain gradient theory the strain energy, W, is assumed to be a function 

W = W(E ij ,d k E ij ) (1) 

of the elastic strain, E^, and the gradient of the elastic strain, dkEij, which is sometimes called 
hyperstrain. The E^ is dimensionless and the d^E^ has the dimension of the reciprocal of 
length. Obviously, no elastic rotation and gradients of it appear in (1). The elastic strain can 
be a gradient of the displacement Ui 

E i:i = d^Uj) = i (diUj + djUi) (2) 

or can have the following form [14-19] 

Eij = d(iUj) + , Eij = Eji (3) 

where denotes the incompatible part of the elastic strain. In addition, the negative in- 
compatible strain may be identified as the plastic strain = The strain (2) is the 
elastic strain in a compatible situation. On the other hand, (3) corresponds to the elastic strain 
for the incompatible case (e.g. dislocations and disclinations) . In the incompatible situation the 
Ui and $ij are not "good" physical quantities because they are discontinuous functions. Thus, 
they cannot be physical state quantities. Because the elastic strain is a physical state quantity, 
it must be a continuous function. It is worth to note that the incompatible strain (3) is similar 
in form to Mindlin's so-called relative strain which contains a micro-strain term. But in our 



2 




Figure 1: Double force stress r^fe acting on a plane with normal rife. 



case this part is identified as the incompatible elastic strain. Thus, we deal with a theory of 
incompatible strain gradient elasticity. 

The most general form of the strain energy for a linear, isotropic, gradient-dependent elastic 
material is given by [3,4] 

W = ^XEuEjj + fiEijEij + c 1 (d j E ij )(d k E ik ) + c 2 (d k E u ) (djE jk ) 

+ c 3 (d k E u ) (d k Ejj) + C4(d k Eij) (d k E l3 ) + c b (d k E l3 ) {d t E 3k ). (4) 

The constants A and /x are the Lame constants and the five c n are the additional constants 
(gradient coefficients) which appear in Mindlin's strain gradient theory [3,4]. Thus, five strain 
gradient terms 1 appear in the isotropic case. 

On the one hand, the Cauchy stress tensor is defined as (Hooke's law) 

dW 

°X? := = X5ijE kk + 2/iEij (5) 

and has the symmetry <n 3 = <Tji. The elastic stress has the dimension of force, Fi, per unit area 
dAj . In addition, the elastic energy is completely symmetric in the strain and stress tensor so 
that the condition for the strain follows 



13 ~ da i3 ~ 2^ V 11 1 + 

with A = 2/iz//(l - 2u). 

On the other hand, the hyperstress tensor is defined as follows (see also [3,4]) 

T v k := \ = c i(SkidiEij + SkjdiEu) + - C2(SkidjEu + S kj diEu + 25i 3 diEi k ) 

0(O k Eij) l 

+ 2c 3 5 t3 d k E u + 2c i d k E lj + c 5 (diE jk + d 3 E ik ) . (7) 

It has the character of double forces per unit area. The dipolar or double force acts through that 
area. In this case r^-fe is a dipolar or double stress tensor. The first index of r^-fe describes the 
orientation of the pair of forces Fi, the second index gives the orientation of the lever arm Axj 
between the forces and the third index denotes the orientation of the normal n k of the plane 
on which the stress acts (see Fig. 1). The double stress gives a contribution in the compatible 
as well as the incompatible case. In strain gradient theory it has the symmetry 

Tijk = T(ij) k = -^{ T i 3 k +Tj ik ), (8) 



1 We want to note that Feynman[20] used five gradient terms of the metric tensor to obtain a linear theory 
of gravity. In that sense the theory of gravity can be considered as a four-dimensional "strain" gradient theory 
in which the metric is the strain tensor and gravitation represents a "metrical elasticity" of space. 
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Fi gure 2: Double force stresses (a) T(y X ^ k and (b) T[y X ^ k acting on a plane with normal n k . 

and possesses 18 components. Thus, the symmetric part T(ij)k arises from double forces without 
moment (see Fig. 2a) . It can be resolved into the dilatational double stress Tuk and the traceless 
symmetric double stress T^j) k — \8ijTu k . No antisymmetric double stresses r^ k , which arise 
from double forces with moment (couple stresses), appear 

T [ij]k = ~ T i %k ) = °' ( 9 ) 

which would correspond to rotation-gradients (see Fig. 2b). In general, a couple stress tensor 
has 9 components. 

For vanishing external body forces, the force equilibrium equation can be derived from the 
principle of virtual work as (variation with respect to the displacement) [3] 

djfaj - d k Tijk) = 0. (10) 

Since we consider an infinitely extended medium, we may neglect additional boundary condi- 
tions. With Eqs. (5) and (7), the force equilibrium (10) reads 

dj^XSijEkk + 2(j,Eij - (ci + c 5 )(didiEij + djdiE H ) - c 2 (d i d j E u + SijdkdiEik) 

-2(c 3 S t3 AE kk + c 4 A£y)) = 0, (11) 

where A denotes the Laplacian. Such an expression (11) is quite formidable and hard to solve. 
Thus, a physical simplification could be useful. 

First possibility of simplification: Let us follow a way in gradient elasticity similar to the 
way which Feynman [20] went in linear gravity. He pointed out that not all five gradient terms 
are necessary. Because the gradient term to C5 can be converted to the C\ gradient term by 
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integration by parts. Therefore, the C5 term may be omitted. We set C5 = 0. In addition, the 
Cauchy stress tensor should fulfill an equilibrium condition of the kind as in elasticity 



djtru = 0. (12) 

Consequently, in (10) it must yield 

dkdjT i: j k = 0, (13) 
and Tijk is self equilibrating. Using (7) and C5 = 0, Eq. (13) reads 

ci (didjdiEij + AdiE H ) + c 2 (diAE u + d l d k diE lk ) + 2(c 3 d t AE kk + c A d 3 AE l3 ) = 0. (14) 
It may be solved by 

didjdiEji(ci + c 2 ) = 0, 
Ad l E il (c 1 + 2c i ) =0, 

Ad i E ll (c 2 + 2c 3 ) =0. (15) 
If we choose a scale such that C4 = 1/2, then the solution of (15) is given by 

I I 

ci = -I, c 2 =t, c 3 = --, c 4 = -. (16) 
The coefficient I has the dimension of a force. Thus, the coefficient I may be chosen 

with only one "gradient-coefficient" or "parameter of nonlocality" 

kT 2 > 0. (18) 

k 2 has the dimension of 1/ [length 2 ]. It is a positive material constant. What we have recovered 
is nothing but the so-called Einstein choice in three dimensions scaled by a factor I (see, 
e.g., [17,21]). So, we can write down the expression 

-jd k T ljk = (incE) y . = -eihi€j mn dkd m Ei n 

= AEij — (djd k E ik + did k E k j) + 5ijd k diE k i + didjE kk — SijAE kk (19) 

which is equivalent to the three dimensional linear Einstein tensor. The linear Einstein tensor 
is the incompatibility of the strain (see [22, 23]). Thus, the three dimensional Einstein tensor is 
proportional to the divergence of the double force stress tensor. The Einstein choice has been 
used in the gauge theory of dislocations by Lazar [17, 18, 24] and by Malyshev [25]. Eq. (19) leads 
to physical solutions for the screw dislocation [17, 18, 25] . The solution of an edge dislocation [25] 
has a modified far field of the stress. In fact, it does not coincide with the far field of the classical 
solution. One way out is to modify the strain energy with additional bend- twist terms [19]. But 
then we are leaving the framework of strain gradient elasticity which is not our aim. Anyway, 
it is interesting to note that the three dimensional Einstein choice which is used in the gauge 
theory of dislocations is contained in the strain energy expression (4) given by Mindlin. 

Second possibility of simplification: If we require that the stress-strain symmetry of the 
elastic energy should also be valid for the gradient terms, and we think this is quite natural, 
we have to use the following choice of the five gradient coefficients 

ci =0, c 2 = 0, c 3 = c 4 = c 5 = 0. (20) 
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Then the strain energy can be rewritten in the following simple form 

W = \ + JL (fan) {dkEij). (21) 

One can see that k has the dimension of an inverse length. Thus, by using the special choice (20), 
we have obtained a simple but still rigorous version of Mindlin's gradient theory which is a 
simple strain gradient as well as stress gradient theory. In particular, the elastic energy (21) is 
symmetric with respect to the strain and the stress and also with respect to the strain gradient 
and the stress gradient 

dW 1 dW 1 

a/a T-i \ o UkVij, -7T77: \ = —n^k^ij- (22) 

d(d k E i:j ) k d(d k (Tij) k 

Then the double stress tensor has the following form 

njk = —dk<y io = — (2fidkEij + \5 io d k Eu). (23) 

Consequently, the double stress (23) is a simple gradient of the Cauchy stress. Thus, it is 
a higher-order stress tensor. Only k~ 2 is a non-standard coefficient of the theory. Recently, 
double stresses similar in their form to (23) have also been used in the mode-I crack problem [26] 
and the mode-Ill crack problem [27,28]. 

Let us mention that a higher order gradient elastic energy with the required symmetry would 
have the following form 

11 1 

w = j a v E ij + ^2 { d k^ij){dkE tj ) + ^ (did k cr i: j) (didkEij) 

+ { d mdid k a lj )(d m did k E i:j ) H , (24) 

where k and k are higher order gradient coefficients. It is interesting to note that the last term 
in Eq. (24) with the third strain gradient corresponds to a force octupole. But in the following 
we only consider the first gradient theory given by (21). 

In addition, we note that Altan and Aifantis [29] have used a similar choice like (20) for the 
study of cracks within their version of strain gradient elasticity. Anyway, they did not use the 
notion of double stress. For that reason, they identified the Cauchy stress with the total stress 
tensor. The price they had to pay was that the singularities are still present in the stresses. In 
order to regularize both the elastic strain and the stress singularities Ru and Aifantis [41] (see 
also [42]) introduced a constitutive relation with gradients of the elastic strain and the stress 
multiplied by two different gradient coefficients. In our framework, which we use in this paper, 
it is not necessary to introduce such a constitutive relation because the Cauchy stress is defined 
as the derivative of the strain energy with respect to the elastic strain and nothing else. On 
the other hand, Gutkin and Aifantis [30] used another choice so that no double stresses appear. 
But, on the other hand, triple stresses which correspond to second gradient strain must appear 
in their approach and, therefore, in the equilibrium equations. But they neglected the triple 
stresses which is not straightforward. The strain gradient elasticity which contains two different 
gradient coefficients proposed by Ru and Aifantis [31] is used by Gutkin and Aifantis [32-34] 
for dislocations and disclinations. 

By substituting (23) into (10), we obtain 

(l - K- 2 A)d j a lJ ; = 0. (25) 

If we define 

(Jjj = (Tij - d k T ljk = X5ijE u + 2fj,Eij - ^(XdtjAEn + 2 / uA£ , ij ), (26) 
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Eq. (25) takes the form 



djOH = 0. (27) 

The stress tensor may be called total stress tensor, which is a kind of "balanced stress" . In 
the gauge theory of defects it is identified as the background stress tensor. One other difference 
with Gutkin and Aifantis' approach [30] is that they consider (26) as a modified constitutive 
relation while we obtained (26) as field equation instead. 

On the other hand, in the incompatible situation we have the additional field Which 
equation corresponds to this field? One should obtain Eq. (26) as a variation of strain energy 
with respect to . But so far this does not work. The way out is to add an additional strain 
energy part to (21), 

W = —(TijEij, (28) 

which in the compatible situation is a null Lagrangian such that it gives no contribution in 
the variation with respect to the displacement. The nontrivial traction boundary problems in 
the variational formulation of field theory can be formulated by means of a null Lagrangian. 
This procedure is well-known in the gauge theory of dislocations [15-19,25]. When the null 
Lagrangian is added to the elastic Lagrangian (strain energy) it does not change the Euler- 

Lagrange equation (force equilibrium) because the associated Euler-Lagrange equation, dj&ij = 
0, must be identically satisfied. Then for the incompatible case, Eq. (26) may be considered as 
the field equation corresponding to /?,.,• . It is interesting that (26) appears in the compatible as 
well as the incompatible situation. Only the interpretation of (26) is different. 

Now we rewrite Eq. (26) and obtain an inhomogeneous Helmholtz equation for every com- 
ponent of the Cauchy stress 

(l-«r 2 A)<7y =aij. (29) 

Since the factor kT 1 has the physical dimension of a length, it defines an internal characteristic 
length in quite a natural way. It is worth to note that Eq. (29) agrees with the field equation 
for the nonlocal stress in Eringen's nonlocal elasticity [35-38] and for the stress in the gradient 
elasticity given by Gutkin and Aifantis [32-34]. If we consider dislocations and disclinations 
which have an axial symmetry, Eq. (29) may be rewritten as a convolution integral 

*«(»■)=/ a(r-r')<^(r')d W (r'), (30) 
Jv 

with the corresponding two-dimensional Green's function 

a(r -r') = ^ K ( Ky /( x - x' f + (y - y'f) • (31) 

Here K n is the modified Bessel function of the second kind and n — 0, 1, . . . denotes the order of 
this function. In comparison with Eringen's nonlocal theory of elasticity the Green function (31) 
may be identified as a nonlocal kernel introduced by Ari and Eringcn [39] . Using the inverse of 
Hooke's law for the stress Oij and Oij , it follows that the elastic strain can be determined from 
the equation 

(l- K - 2 A)E ij =E ij , (32) 

o 

where E^ is the classical strain tensor. Because the strain tensor fulfills an inhomogeneous 
Helmholtz equation, we may rewrite (32) as a nonlocal relation for the strain 

E ij (r)= I a {r-r') E l3 {r')Av{r'). (33) 
Jv 
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Of course, such a relation (33) does not appear in Eringen's theory [38] of nonlocal elasticity. In 
his theory, the displacement and the elastic strain are the same as in classical elasticity. Using 
Hooke's law, we can combine (29) and (32) to a gradient like relation 

(1 - k- 2 A)^ = (1 - K- 2 k)(\5 t] E kk + 2jiE t] ). (34) 

It is interesting to note that Eq. (34) has the same form as the gradient constitutive relation 
given in [31] if their two different gradient coefficients are the same gradient coefficients. For 
modified solutions of the stress and elastic strain fields we require that the far field of them 
should agree with the classical expressions and they should be free from the classical singularities 
at the defect line. Therefore, we are looking for nonsingular solutions for both the stress and the 
elastic strain. It depends on the taste of the reader to consider alternatively such constraints 
as physically motivated boundary conditions. 

Using the decomposition (3), we obtain the coupled partial differential equation 

(1 - «- 2 A) [d {iUj) = d (l u 3) ^. )5 (35) 

where Ui denotes the displacement field and is the plastic distortion in classical defect theory 
(see, e.g., [14]). Thus, if the following equations are fulfilled 

(l- K - 2 A)/3y =kj, (36) 

(l-«- 2 A)/3£=fe (37) 

the equation for the displacement field 2 , 

(l-/s- 2 A)« i = u i) (38) 

is valid for the incompatible case. Thus, for defects (dislocations, disclinations) the inhomogc- 
neous parts of Eqs. (37) and (38) are fields with discontinuities (jumps). We note that Eq. (38) 
was used by Gutkin and Aifantis [32,43,44] in order to calculate the displacement fields for 
screw and edge dislocations. 

In nonlocal elasticity Eringen [35-38] found the two-dimensional kernel (31) by giving the 
best match with the Born-Karman model of the atomic lattice dynamics and the phonon dis- 
persion curves. The length, k _1 , may be selected to be proportional to the lattice parameter a 
for a single crystal, i.e. 

= e a, (39) 

where eo is a non-dimensional constant [36] . Obviously, for eo = we recover classical elasticity. 
Eringen [35,36] used the value of eo = 0.39 in nonlocal elasticity. 

We notice that a negative gradient coefficient kT 2 < would change the character of the 
Hclmholtz equations (29) and (32) and the corresponding solutions. Let us emphasize that the 
solutions which we consider in the following sections are valid for a positive gradient coefficient. 
In addition, non-negative definiteness of the strain energy (21), W > 0, requires the following 
conditions for the material constants and the gradient coefficient (sec, e.g., [45]) 

3A + 2^>0, \x > 0, «T 2 > 0. (40) 

3 Dislocations 

In this section we consider straight dislocations whose line coincides with the z-axis of a Carte- 
sian coordinate system in an infinitely extended medium. 

2 If /3fj = (compatible distortion), the inhomogeneous Helmholtz equation, which was already proposed by 
Aifantis [40], Ru and Aifantis [41], (l - k~ 2 A) Ui = Ui is obtained without further assumptions. 
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3.1 Screw dislocation: 6 = (0, 0, b z ) 

We start with the simplest case, the anti-plane strain which corresponds to a screw dislocation. 
We make an ansatz for the total stress and for the Cauchy stress which has the form as 



an = 



( -d y F \ 

d x F 
\ -d y F d x °F / 



-d y F 



<nj = d x F . (41) 



-d y F d x F 



We choose for F the well-known stress function of elastic torsion, sometimes called Prandtl's 
stress function. It is given by (see, e.g., [23]) 

F=^lnr, (42) 

with r = \J x 1 + y 2 . Substituting (41) and (42) into (29) we obtain for the stress function F 
the following inhomogencous Hclmholtz equation 

(l- K - 2 A)F=^ lnr. (43) 

The nonsingular solution of (43) is (see, e.g., [16, 18,37]) 

F=^{lnr + Jf (/sr)}, (44) 

which represents a stress function for a nonsingular screw dislocation. In the far field, the stress 
function (44) agrees with Prandtl's stress function and for small r it cancels the logarithmic 
singularity. Consequently, the elastic stress is given by 



the corresponding field of elastic strains reads 



E zx — 

47r r 



^{l-Krtfi(/tr)}, E zy = ^-^{l- K rK 1 ( K r)}. (46) 



The appearance of the modified Bessel function in (45) and (46) leads to the elimination of 
classical singularity ~ at the dislocation line (see Eq. (A. 2) in Appendix A). The stress <r zy 
has its extreme value \a zy (x, 0)| ~ 0.399k-^ at |a:| ~ 1.114k -1 , whereas the stress a zx has its 
extreme value \a zx (0,y)\ ~ 0.399k^ at \y\ ~ 1.114k _1 . The strain E zy has its extreme value 
\E zy (x,0)\ ~ 0.399k|^ at \x\ ~ 1.114k- 1 and E zx has its extreme value \E zx {Q,y)\ ~ 0.399k|^ 
at \y\ ~ 1.114k _1 . In addition, the stresses and strains are zero at r = 0. The stress is 
plotted in Fig. 3. It is worth to note that the stress (45) agrees with the nonlocal stress given 
by Eringen [35-38], with Edelen's expression [16] and with the stress given by Gutkin and 
Aifantis [32-34]. 

Let us now give some remarks on the double stresses of the screw dislocation. The non- 
vanishing components are given by (and the components due to the symmetry T^jjk) 

T zyx = d xx F, r zxy — — — d yy F, T zyy = ~t zxx = — d xy F. (47) 

r\j r\i rv 

They have a similar form like the elastic bend-twist tensor of a screw dislocation given in [46] . 
This means that the double stresses are singular at r = 0. Because the double stress is a simple 
gradient of the Cauchy stress which is not singular in our case, it is less singular as a gradient 
of the stress calculated in classical elasticity. 
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(a) 



(b) 



Figure 3: Stress of a screw dislocation near the dislocation line: (a) a xz and (b) a yz are given 
in units of \ib z nj\l^\ . 



3.2 Edge dislocation: b=(b x , 0,0) 

In the case of plane strain, we may make the following stress function ansatz 

/ o o \ 

1 d$ y f o -dlyf > 

£«= -dlj dlj , <m = 

V uAf J 



dU 'dlj 

-d 2 f d 2 f 

^xyJ ^xxJ 

vAf 



(48) 



Obviously, it yields a zz — v{o xx + o yy ). For an edge dislocation with Burgers vector b x , we use 
the corresponding Airy's stress function [23] 

^b x 



/ = -; 



• j/lnr. 



(49) 



2tt(1 - i/) 

If we substitute (48) and (49) into (29), we get the inhomogeneous Helmholtz equation for the 
stress function / 

(J-b x 



(l-K- 2 A)/ 



y In r. 



(50) 



2tt(1 - v) 

The nonsingular solution for the modified stress function of a straight edge dislocation is given 
by [19,46] 



/ = 



[ib x 



y< lnr + 



2tt(1 - v) 

By means of Eqs. (48) and (51), the elastic stress is given as [19] 
lib 



(51) 



®xx — 



2tt(1 



^ 2tt(1 - v) 



r^) ^{V + 3x2 ) + ^2 - 3a;2 ) - 2y 2 K rifi( K r) - 2(y 2 - 3a: 2 )X 2 M J, 

K, T J 

2 y 2 ) - ^(x 2 - 3y 2 ) - 2y 2 K rK 1 ( K r) + 2(x 2 - 3y 2 )^ 2 ( K r)|, 



' ra 2tt(1 -is)r^ [V X 



x 



a zz = —^-^{l- K rK l{K r)\. 
7r(l — v) r z I ) 



(52) 
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Figure 4: Stress of an edge dislocation: (a) <j XXl (b) a xy , (c) a yy are given in units of 
ldb x K/{2Tr(l — v)] and (d) a zz is given in units of fib x VK/[n(l — v)] . 



The stress (52) is zero at r = 0. In fact, the "classical" singularities are eliminated due to 
the behaviour of the Bessel functions at the dislocation line (see Appendix A). The stress (52) 
has the following extreme values: \a xx (0,y)\ ~ 0.546k 2 i!(i-v) at If I — 0.996k- 1 , \<J yy (0,y)\ ~ 
0.260^2^^ at \y\ ~ 1.494k- 1 , ^(a;, 0)| ~ 0.260*2^^ at \x\ ~ 1.494k- 1 , and |a„(0,y)| ~ 
0.399/t at |y| ~ 1.114k _1 . The stress is plotted in Fig. 4. The corresponding trace of the 
stress tensor a k k = &xx + &yy + &zz produced by the edge dislocation is 

a kk = -^f^^{l-nrK l{ nr)\. (53) 
7T(1 — v) r z L J 
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Using the inverse of Hooke's law, we find for the elastic strain of this edge dislocation 

^ b x y f,, „ , 2x 2 4 , , „ 

47r(l — v) r z \ r A K-^r 4 ' 



- 2 (j^ - nrK^nr) - ^{y 2 - 3x 2 )K 2 {nr) 



(54) 



x 2 



- 2 

,2 



^ Kr^i(Kr) + ^(y 2 - 3x 2 )K 2 ( K r) J, 



^ = -|^^{l-^M}- (55) 



The strain (54) has the extreme values {v = 0.3): \E XX (Q, y)\ ~ 0-308K; 3¥ ^ r ^ at \y\ ~ 0.922k" 1 , 
|£ w (0,y)| ~ O.OIOk^J^ at |y| ~ 0.218k- 1 , |£? ra (0,y)| ~ 0.054^^-^ at \y\ ~ 4.130k" 1 , 
and |£J xl/ (a:,0)| ~ 0.260K^j-^y at a; ~ 1.494k" 1 . In addition, it is interesting to note that 
E yy (0,y) is much smaller than E xx (0,y) within the core region (see also [32]). The dilatation 
Ekk reads 

b x (l-2v) y 
2tt(1 - v) r 2 

The strain and the dilatation are zero at the dislocation line. 

It is interesting to note that the stress (52) and the strain (54) of a dislocation with Burgers 
vector b x agree with the expressions calculated by Gutkin and Aifantis [32-34] by using the 
technique of Fourier transformation. 

3.3 Edge dislocation: 6=(0, &j,,0) 

Let us now complete the case of dislocations with the edge dislocation with Burgers vector b y . 
Again it is a plane strain state such that we can use the stress function ansatz (48). But now 
Airy's stress function reads 

Substituting (48) and (56) into (29) we find 

flby 



(^"-^'-sfe* 1 "- (57 » 

Its solution is given by 



(58) 



Then we find for the elastic stress 



> xy 



2tt(1 - 


v) r 4 


flby 


X 


2tt(1 - 


v) r 4 


liby 


y 


2tt(1 - 


v) r 4 


flbyV 


X f 



(x 2 + 3y 2 ) + ^{x 2 - 3y 2 ) - 2x 2 KrK 1 ( K r) -2{x 2 - 3y 2 )K 2 ( K r) 
{x 2 - y 2 ) - ^j(3^ 2 - y 2 ) + 2x 2 K rK 1 ( K r) + 2(3x 2 - y 2 )K 2 (Kr) J, 
r(1 _ i/)rar -*rJf 1 («r)}, (59) 
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and for the trace of the stress tensor 

"--'H^^ 1 (») 

The stress (59) has following extreme values: |cr yy (x,0)| ~ 0.546k 2 i!(i V -v) a ^ M — 0-996k _1 , 

|a xx (a;,0)| ~ 0.260 K 2i^y at \x\ ~ 1.494k- 1 , |<7 xy (0,y)| ~ 0.260 K ^g^j at |y| ~ 1.494k- 1 , 

and |cr^(x,0)| ~ 0.399k^^j at |a;| ~ 1.114k" 1 . 
The elastic strain is given by 



4tt(1 



r 2 K 2 r 4 



~ 2 (^2 - ") rf i( Kr ) + ^(^ - 3y 2 )K 2 (Kr)J, (61) 



E yy = ITT — T 4( a - 2 -) + ¥ + ^4 (- 2 - 3 y 2 ) 

47r(l — f) H ^ r A K^r 4 v 



a: 2 \ 2 1 

-2 - i/J Kr^Kr) - — (a; 2 - 3y 2 )if 2 (KrU, 



1-Jtai — 



and the dilatation reads 



^-sn^pf 1 -"*■<">}• (62) 

The strain (61) has the extreme values (1/ = 0.3): \E yy (x 1 0)\ ~ 0.308k 47r ^_ ;y) at |x| ~ 0.922k- 1 , 
|£ xx (a;,0)| ~ O.OlOK^p^y at |z| ~ 0.218k- 1 , \E xx (x,0)\ ~ 0.054/e 3;F ^_ y at |z| ~ 4.130k- 1 , 
and \E xy (0, y)\ ~ 0.260kj^j^j at ~ 1.494k -1 . Again, it is interesting to note that E xx (x, 0) 
is much smaller than E yy (x,0) within the core region. In addition, the stress and strain holds 
are zero at r = 0. 

In principle, we could calculate the double stresses of edge dislocations. But we do not want 
to do this in detail. Using Eqs. (52) and (59) we would obtain expressions which are similar in 
form to the double stresses of the screw dislocation. Again, the double stresses are still singular 
at the dislocation line. The double stresses 3 of an edge dislocation are given in terms of the 
stress function / as derivatives of the third order according to: 

T yyx — — 2 d xxx f, T yyy = ~ T xyx — — d xxy f 

K K 

T XX y = ^2 dyyyf , = ^T^y = ~^ ^y^ 

- V{r xxx -\~ Tyy X ^ T Z zy V{jxxy ~t~ ^yyy) . (^3) 



They are singular at r = 0. 



4 Disclinations 

In this section we consider straight disclinations in an infinitely extended medium. The discli- 
nation line coincides with the z-axis of a Cartesian coordinate system. We are using deWit's 
expressions [14] for the classical stress and strain fields (see also [48]). 

3 In the meantime [47], we have calculated the double and triple stresses of screw and edge dislocations in 
second strain gradient elasticity. The double stresses can be found there in the limit from second to first strain 
gradient elasticity. 
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4.1 Wedge disclination: CI = (0, 0, Q z ) 

As in the case of edge dislocations we may use the stress function ansatz (48) for the wedge 
disclination as well. It is obvious because the wedge disclination corresponds also to plane 
strain. Using the stress function of a "classical" wedge disclination, 

/ifl z 2 f i 1 — 4^ 



/= 4^r r r -w^)\> (64) 

and (48), we obtain the following equation which gives the solution of a "modified" stress 
function of a wedge disclination 

f 1 - «~ 2 a) f = ~T~T~~ — \ r 2 (inr ^^-) . (65) 
V ) J 4?r(l — i/) V 2{l-2v)J v ' 

Consequently, the "modified" stress function of a wedge disclination is given by (see also [49] ) 



Substituting (66) into (48), we find for the stresses of a wedge disclination 

f 1 , y 2 , v , is , s , i x2 - v 2 ) 



2tt(1 - v) 



{ Inr + V y 2 + ^ + *o(«r) + i^-fi (2 - .W 2 (.r) 

2^) { 4 + + - ^ (2 - *W 2 (.r) 

( 2 - K 2 r 2 A 2 ( K r)) J, 



^ _ 2tt(1 - 1/) r 2 i ' 



tt(1 — 1/) I 2(1 -2i/) 

and for the trace of the stress tensor 



a kk = !±L-> ( Inr + ^-^ + ^oH[ (68) 

The stress is plotted in Fig. 5. Using Eqs. (A.l) and (A. 3), we obtain for the stress at r = 

^(0) = a yy (0) = - 7 - In £ + I}, ^(0) = 0, 

*.,(0) = ^-a fcfc (0) = -f^ {-V - 7 - In f j. 
1 + v 7r(l — v) \ \ — 2v 2 J 

Consequently, the stress is finite at the disclination line in contrast to the unphysical stress 
singularity in "classical" disclination theory. The elastic strain is easily calculated as 

E xx = ^^{(1 - 2„)(]nr + *o(«r)) + ^ + (2 - «V* a («r)) }, 

^ = {(1 - 2.) (Inr + A oM ) + £ - (2 - «W a («r)) }, 

{l--^( 2 - K VA' 2 M)}. (69) 



xv ~ 4tt(1 - 1/) r 2 
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Figure 5: Stress of a wedge disclination: (a) a xx , (b) a xy , (c) a yy are given in units of 
/xf2 z /[27r(l — v)\ and (d) <j zz is given in units of /xf2 z i//[7r(l — ^)] . 



and for the dilatation we obtain 

Ekk = 2AI- v) { (1 ~ 2V) ( ln r + K ° {Kr) ^ + (?0) 
Again using Eqs. (A.l) and (A. 3) the strain reads at r = 

E xx {0) = E yy (0) = l -E kk (Q) = {(1 - 2i/)( 7 + ln£) - 1|, £ xy (0) = 0. 

Now we calculate the non-vanishing components of the double stress for the wedge disclina- 
tion. They are given in terms of the stress function as follows 

T yyX = — 2 d xxx f, T XXy = — d yyy f, Tyyy = —T XyX = — d yxx f, 

1% rv rv 

T X xx — ~ T xyy — ~ 9 xyy f, T zzx = — 9 X A/, T zzt/ = ~ 9 y Af. (71) 

fv rC rv 
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So we obtain for the double stresses of the wedge disclination 

4 /0 „ c 



Tyyx — 



' xxy 



+ 3y 2 ) 



'yyy 



7~zzx — 



2 - 3y 2 ) - 2x 2 KrK 1 (nr) -2(x 2 - 3y 2 )K 2 (Kr)^, 
^^{(V 2 + 3- 2 ) + ^(V 2 - 3- 2 ) - 2 y 2 ^ l( Kr) - 2( y 2 - 3.* 2 )X 2 m}, 
^ |V - x 2 ) --L^(y 2 ~ 3x 2 ) - 2x 2 nrK 1 {nr) + 2(y 2 - 3x 2 )K 2 ( K r)^, 





X 


2ir(\ - v\k 2 






;<y 


2tt(1 - v)k 2 






;<y 


2tt(1 - v)k 2 


r 4 




X 


2tt(1 - f)*; 2 


— 




1/ 


7r(l — i^)k 2 r 


2 I 




' r 


7r(l — ^)k 2 r 2 I 



(72) 



This double stress tensor has some interesting features. First, it is not singular. Second, it has 
a similar form like the stress field of an edge dislocation (compare with Eqs. (52) and (59)). 
Thus, the components are zero at r — and have extrcmum values near the disclination line. 

4.2 Twist disclination: fi = (Q x ,0,0) 

In the case of twist disclination the problem, which we want to consider, is more complicated 
than those of dislocations and of the wedge disclination. The reason is that the situation is no 
longer a proper two-dimensional problem. In the case of twist disclination the three-dimensional 
space may be considered as a product of the two-dimensional xy-plane and the independent 
one-dimensional 2-line [14]. Thus, the z-axis plays a peculiar role. 

First, we make an ansatz which fulfills the stress equilibrium. It is given by 





I 




d 2 f 




~d 2 y f 


-d v F + d z °g 


\ 




Gij = 






-d 2 f 


d 2 J 


d x °F 






\ 


-d y F + d z °g 


d x °F 


O 

P 


J 






d 2 f 

yyJ 


-d 2 f 

^xyJ 


~d y F + O z g 




(Tij 


! 


-d 2 f 


d 2 f 


d x F 


) 






-d y F + 8 z g 


d x F 


P 





with the relations 



P = fA/ = -d x g, p = vAf = -d x g, 



(73) 



(74) 



which follow from did z i = and diO Z i = 0. One can see the special role of z in the ansatz (73). 
The ansatz (73) is not only an addition of the anti-plane (41) and plane strain (48) situation 
because an additional stress function g or g enters the ansatz. The following "classical" stress 
functions, 



/ = 

O 

F 



2tt(1 - v 
[£l x 



: In 7 



2tt(1 - v) 
y£L x v 



ylnr, 



7T(1 - V) 



zln r, 



(75) 



1G 



may be used to reproduce de Wit's expressions for the stress and strain fields of the twist discli- 
nation. Substituting (75) and (73) into (29), the following Hclmholtz equations to determine 
the "modified" stress functions follow 



tt(1 - v) 

The solutions of the modified stress functions are given by 



2tt(1 - v) 



2tt(1 - j/) 



r(l-f) 

By means of (73) and (77) we are able to calculate the stress 



fe) 5^ - ^ - k4^ - 3y2 ) - V«r^i(«r) + 2(x 2 - 3 y 2 )^ 2 (Kr)}, 



2tt( 
2tt(T 



- ^|(x 2 + 3y 2 ) + ^(z 2 - 3y 2 ) - 2x^rK 1 (nr) 2(x 2 3y 2 )^ 2 ( K r)|, 

(x 2 - y 2 ) - -^(3.x 2 - y 2 ) + 2x 2 K rK 1 ( K r) + 2(3x 2 - y 2 )K 2 (nr) ), 
K, T J 



2tt(1 - v) r 4 

^ X f (l - 2„)(lnr + *>(«r)) + £ + ( 2 - K 2 r 2 K 2 ( K r)) }, 



2tt(1 - 

'--sjr^?! 1 -^^-^"")}- (78 » 

The trace of the stress tensor reads in this case 

ni = _^A^ {l _ KrKi{Kr) y (79) 

The stress (78) has its extreme values in the xy-planc: \a yy (x, 0)| ~ 0.546k 2 ^i- v ) at \ x \ — 
0.996k" 1 , \a xx {x,0)\ ~ 0.260^2^^ at |z| ~ 1.494k- 1 , \a xy {0,y)\ ~ O^K^^y at |y| ~ 
1.494k- 1 , and |t7 zz (a;,0)| ~ 0.399k ^^ at |x| ~ 1.114k -1 . The stress a zx has at r = the 
value: a z;c (0) ~ 5^j[(l - 2i/)( 7 + In f ) - 0.5] and with v = 0.3: ^(0) ~ ^^[0.4 In k - 
0.546]. 
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The corresponding elastic strain is given by 

Q x zx ( , . 2y 2 4 , 9N 

,2 



^ - -4^) ^ ~ ^ + ^ + ^ ^ } 



^2 - i/J nrKi(nr) - -^(x 2 - 3y 2 )K 2 (nr) 

E*y = f {i - |r + ^i(3- 2 -y 2 ) - |r-M-) - ^(3, 2 - y 2 )* 2 (.r)}, 



(80) 



The dilatation reads 



The strain (80) has the extreme values (v = 0.3): \E yy (x 1 0)\ ~ 0-308 ^"^ at |x| ~ 0.922k" 1 , 

at |z| ~ 0.218k- 1 , \E xx (x,0)\ ~ 0.054^ 
and |-E X!/ (0, y)| ~ 0.260k 4 9^ , at |y| ~ 1.494k- 1 . Again, it is interesting to note that E xx (x, 0) 



\E xx (x,Q)\ ~ O.OlOK^^y at |z| ~ 0.218k- 1 , |£ xx (z,0)| ~ 0.054k s ^- j at |z| ~ 4.130k" 

t(1 

is much smaller than E yy (x,0) within the core region. The strain E zx has at r = the value: 
E zy (0) ~ 47r ^s [(1 - 2i/)(7 + lnf ) - 0.5]. The dilatation has its extremum \E kk (x,0)\ ~ 

0.399k%^§^ at |z| ~ 1.114k" 1 . 

It is interesting to note that there is a relation between the stress and strain fields of twist 
disclinations and edge disclinations. If one replaces £l x z by —b y in the components of the stress 
<?xx ~ &zz in (78) and in the strain E xx - E xy in (80), the stress (59) and the strain (61) are 
obtained. 



4.3 Twist disclination: Q, = (0,^,0) 

For the twist disclination with Frank vector £l y we make a similar ansatz like the ansatz (73) 
of the twist disclination with Frank vector Cl x . It is given by 

-d 2 y f -d y F \ 

dU 8 x F + d z g , 
d x F + d z g p J 

(82) 

with the relations 

P = vAf = -d v g, p = uAf = -d v g, (83) 

which follow again from di<7 Z i = and di<j Z i = 0. The only one difference between (73) and 
(82) is the position of the stress function g. In the present case, the "classical" stress functions 



( 



d 2 f 



yy-> 
-d 2 f 



dl.f 



d 2 f 

xyJ ^xxJ 

V -dyF d x F + d z g 



-dyF 

d x F + 8 z g 

o 

V 



&ii = 



8 2 f 

w yyJ 
-dyF 
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are given by 

° fifty 

2tt(1 - i/) 

g= y zlnr. 84 

7r(l - V) 

These three stress functions reproduce the stress and strain fields given by deWit [14]. If we 
substitute (82) and (84) into (29), we get the following inhomogencous Hclmholtz equations 



fifty 



x In r, 



2tt(1 - j/) 



'(l-x) 

The solutions of the "modified" stress functions reads now (see also [50] ) 



f Ji , 2 

J = — - — j- — - — r zy< mr H — — 



2.(1-,) ~n ^-(l-rc^r))}, 
7r(l — f) L > 



r(l-iz) 

So we find for the elastic stress of this disclination [50] 

fitly zy 

° xx = ~2^r^) t~ 

-2^) 3{ (2/2 - - ^ (?/2 - 32;2) - 2 * 2 ^M + 2(y 2 - 3* 2 )X 2 ( K r)} 
f^T) ^{^ ^ -^ 2 - - 2 ^iM + 2(, 2 - 3 y 2 )^ 2 ( K r)}, 



^||(y 2 + 3a; 2 ) + ^(y 2 - 3x 2 ) - 2y 2 K rK 1 ( K r) - 2(y 2 - 3x 2 )X 2 ( K r) j, 



>xy 



&zz = 



2tt( 



'zy 



2tt(1 

and its trace is given by 



- {(1 - 2,)(lnr + ff («r)) + £ - (2 - K 2 r 2 K 2 (nr)) }, (87) 



^ = -^#^^{l-^ l(K r)}. (88) 
7r(l — ,) L J 



It can be seen that the stresses have the following extreme values in the zy-plane: \cr xx (0, y)\ ~ 
°- 546K 27(fey at Iwl - 0.996k" 1 , \a yy {0,y)\ ~ 0.260^2^^ at l»l - 1-494k"\ |a^(x,0)| ~ 
°- 260K 2^fe) at M - 1.494k- 1 , \a zz (0,y)\ ~ 0.399 K ^^ at |y| ~ 1.114k" 1 and |a fcfc (0,y)| ~ 
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(e) (f) 

Figure 6: Stress of a twist disclination: (a) a xx , (b) a xy , (c) a yy are given in units of 
fiQyKz /[2n(l — v)], (d) a zz is given in units of fxClyUKZ — u)\, (e) a zx and (f) a zy are 
given in units of fj,Cl y /[2tt(1 — v)] . 
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0.399k ^%xz7T^ a * 1 2/ 1 — 1-H4k l . The stresses a XXl a yy and cr^ are modified near the 
disclination core (0 < r < 12k- 1 ). The stress a zz and the trace Okk are modified in the 
region: < r < 6k . Far from the disclination line (r ^> 12k -1 ) the modified and the 
classical solutions of the stress of a twist disclination coincide. In addition, it can be seen that 
at z = the stresses a xx - a zz are zero. The stress a zy has at r = the value: cr zy (0) — 
2n£- v ) 10 ~ 2^)(7 + In §) - 0.5] and with v = 0.3: a zy (0) ~ g^rr^ [0.4 In k- 0.546]. The stress 
is plotted in Fig. 6. 

For the elastic strain we obtain 



47r(l - v) r 2 



y 2 \ 2 1 
v\ KrK\{nr) ^(y 2 — 3x 2 )K2(Kr) >, 



„ Q v zy f ,„ „ . 2a; 2 4 - 9 

f» 47r(l - i/) r 2 { y ' r 2 n 2 r 4y 



2\ 



x 2 



i?j KrKi{nr) + -^(y 2 - 3x 2 )K 2 (nr)^, 



E*v = f {l - - ^ - V) - % nrK\ (kt) + ^(x 2 - 3y 2 )K 2 ( K r) 



l K 2 r 2 ( 



(89) 



The dilatation is 



Again, if one replaces fl y z by 6 K in the components of the stress a xx - a zz in (87) and in the 
strain E xx - E xy in (89), the stress (52) and the strain (54) is reproduced (see also the discussion 
in [50]). The components of the strain tensor have in the xy-plane the following extreme values 
{y = 0.3): \E xx (0,y)\ ~ 0.308^4^^ at \y\ ~ 0.922k- 1 , \E yy (0,y)\ ~ O.OlOK^^y at 

\y\ ~ 0.218k" 1 , |£ TO (0, y)| ~ 0.054*4^^ at \y\ ~ 4.130k- 1 , and |£ xy (z, 0)| ~ 0.260*5$^ 
at a; ~ 1.494k- 1 . It is interesting to note that E yy (0, y) is much smaller than E xx (0, y) within 
the core region. The strain E zy has at r = the value: E zy (0) ~ 4^3^ [(1 — 2^)(7+ln |) — 0.5]. 

The dilatation has its extremum \Ekk(0,y)\ — 0.399k n ^(i"^j Z at \y\ — 1.114k -1 . 

It is interesting to note that the solutions of stress and strain fields of disclinations given in 
this section agree with the expressions earlier obtained by Gutkin and Aifantis [33, 34] by using 
the technique of Fourier transformation. 

Now some remarks on the double stresses of twist disclinations are in order. Only the 
components of double stresses of twist disclinations t xxz , r yyz , T xyz , t zzz , t zxx , T zxy , r zyx and 
r zyy are nonsingular. Again, they are similar in the form like the stresses of edge dislocations. 
The other components t xxx , T yyxi T xyx , t ZZXi T xxy , T yyy , T xyy and r zzy are singular at r = like 
the double stresses of edge dislocations. Moreover these components are zero at z — 0. 

5 Conclusions 

We investigated two special versions of first gradient elasticity. One special version of Mindlin's 
gradient theory has been used in the consideration of dislocations and disclinations. This theory 



21 



is a strain gradient as well as a stress gradient theory. We have applied this theory to all three 
types of dislocations and disclinations. Using the stress function method, we found modified 
stress functions for dislocations and disclinations. These stress functions are modifications of the 
classical ones (e.g. Prandtl's and Airy's stress functions). All modified stress functions satisfy 
two-dimensional inhomogeneous Helmholtz equations in which the classical stress functions are 
the inhomogeneous parts. Using these modified stress functions, exact analytical solutions for 
the elastic stress (Cauchy stress) and the elastic strain fields of all six types of Volterra defects 
have been found. These stress and strain fields fulfill inhomogeneous Helmholtz equations 
in which the inhomogeneous parts are given by the classical singular stress and strain fields, 
respectively. The main feature of these solutions is that the unphysical singularities at the 
defect line are eliminated. Thus, the improved stress and strain fields have no singularities 
in the core region unlike the classical solutions of defects in elasticity which are singular in 
this region. The stresses and strains are either zero or have a finite maximum value at the 
defect line. The maximum value of stresses may serve as a measure of the critical stress level 
when fracture and failure may occur. In addition, if we equate the maximum shear stress to the 
cohesive shear stress, one can obtain conditions to produce a dislocation or disclination of single 
atomic distance. In gradient elasticity the maximum value of stresses depend on the gradient 
coefficient k. Thus, one could test the value of k if one uses the theoretical shear stresses based 
on lattice dynamics calculations. 

The gradient theory considered in this paper contains double stresses (hyperstresses) which 
are simple gradients of the Cauchy stress. In the case of dislocations these double stresses are 
still singular at r = 0. Only for the wedge disclination the double stresses are nonsingular. For 
the twist disclinations some components of the double stress tensor are nonsingular and the 
other components have singularities at r = 0. 

Finally, we notice that the stress fields of dislocations and disclinations are also solutions in 
Eringen's nonlocal elasticity by using the nonlocal kernel (31) and the condition (12). In fact, 
the nonsingular stresses correspond to the nonlocal ones. Eventually, the nonsingular stresses 
of dislocations and disclinations can be calculated by means of the convolution of the singular 
classical stresses with a nonlocal kernel which is a kind of a distribution function. It is the 
convolution with a suitable kernel that provides smoothing of the "classical" stress singularities 
and produces nonsingular stresses which are the main features in nonlocal elasticity. 
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A Expansion of Bessel functions 

In this appendix, we give the expansion of modified Bessel functions which we need to study 
the near field of nonsingular stresses and strains. The expansion is given by (see, e.g., [51]) 



K (nr) 



, nr 
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(A.2) 
(A.3) 
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where 7 = 0.57721566... is Euler's constant. The first terms in the expansions (A.1)-(A.3) 
eliminate the classical singularities of the stresses and strains. This elimination of singularities 
is a kind of "regularization" of stresses and strains. The other terms are zero at r — 0. Thus, 
the second terms in (A.1)-(A.3) describe the behaviour near the defect line in the first order of 
the expansion. 
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